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Abstract 

We report a preliminary measurement of the slope parameter a for the rj Stt" decay carried out with KLOE 
at DA<1?NE; where a is the parameter describing the energy dependence of the square of the matrix element for 
this decay. By fitting the event density in the Dalitz plot with a collected statistic of 420 pb~^ we determine a = 
—0.027 ± 0.004 (stat) Iq ooe (sj/si)- This result is consistent with current chiral perturbation theory calculations 
within the unitary approach. 



1. Introduction 



kinematical variable: 



With a branching ratio of 32% the rj — > Stt^ de- 
cay is a major decay mode of the r] despite the fact 
that is a G-parity forbidden transition. Neglect- 
ing a small electromagnetic transition (Suther- 
land's theorem [T]) this decay is due almost ex- 
clusively to the isospin breaking part of QCD: 



£/ = -- (m„ - md) {uu - dd) 



(1) 



and provides a nice way to determine the up-down 
quark mass difference. 

Many theoretical issues [2], [3] refer to the study 
of the charged and neutral decay mode of ry ^ Stt, 
and in particular to the Dalitz plot parameters of 
these decays. 

The Dalitz plot distribution of the 77 — ^ Btt^ de- 
cay is conventionally described in terms of one 
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It. 



ZEi — rriri 
Tin — 3m^o 



Pmax 



(2) 



where Ei denotes the energy of the i-th pion in 
the r] rest frame and p is the distance from a point 
on the Dalitz plot to its center, pmax is the max- 
imum value of p. For the decays into three iden- 
tical particles, it is possible to use a symmetrical 
Dalitz plot where the event density is described 
by a single quadratic slope parameter a, which 
represents the difference from pure phase space: 



I A^3^« [z) P ~ 1 + 2az. 



(3) 



The lowest order predictions of Chiral Pertuba- 
tion Theory quote a zero value for a. The event 
density in the Dalitz plot must be uniform. A 
non zero value is instead expected by dispersive 
calculation [3] where the effect of tt — tt rescat- 
tering is included. Lately a theoretical result [S], 
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obtained in the chiral unitary approach based 
on the Bethe-Salpter equation, provide for a 
the value —0.031. There are three previous ex- 
perimental determination of a: the GAMS 2000 
group i] quoted a = -0.022 ± 0.023 based on 
5 X 10^ events; the Crystall Barrel Collaboration 
obtained a = —0.052 ± 0.020 from a sample of 
10 X 10* events and finally the Crystal Ball [8] 
result, a = -0.031 ± 0.004 based on 10^ events. 
This was the first statistically significant mea- 
surement of a which agree very well with the 
most recent theoretical result. In this paper 
we report on a new precise measurement of the 
Dalitz plot parameter for the 77 —> Stt^ decay. 



2. DA$NE and KLOE 

The DA$NE e+e^ collider operates at a to- 
tal energy W = 1020 MeV, the mass of the 
0(lO2O)-meson. Approximately 3 x 10^ 0-mesons 
are produced for each pb~^ of collected lumi- 
nosity. Since 2001, KLOE has integrated a to- 
tal luminosity of about 2.5 fb^^. Results pre- 
sented in this paper are based on data collected 
on 2001 — 2002 only and correspond to about 
420 pb^^. The KLOE detector consists of a large 
cylindrical drift chamber, DC, surrounded by a 
lead/scintillating- fiber electromagnetic calorime- 
ter, EMC. The drift chamber [9], is 4 m in di- 
ameter and 3.3 m long. The momentum reso- 
lution is a{pT)/pT ~ 0.4%. Two track vertices 
are reconstructed with a spatial resolution of ^ 
3 mm. The calorimeter |10j, composed of a bar- 
rel and two endcaps, covers 98% of the solid an- 
gle. Energy and time resolution are a{E)/E — 
5.7%/^/E[GeY] and a{t) = 57ps/^/E[GeY] © 
100 ps. A superconducting coil around the de- 
tector provides a 0.52 T magnetic field. The 
KLOE trigger [11], uses calorimeter and drift 
chamber information. For the present analysis 
only the electromagnetic calorimeter (EMC) sig- 
nals have been used. Two local energy deposits 
above threshold, Eth > 50 MeV for the barrel and 
-Eth > 150 MeV for the endcaps, are required. 



3. Dalitz plot of ry ^ tt" 7r° vr" decay 

At KLOE the 77 meson is produced in the pro- 
cess (j) —> r]"f where the recoil photon (E^ = 
363 MeV) is monochromatic and easily selected. 
Thus to select the final state we require to have 
seven prompt clusters in the event. After apply- 
ing a kinematic fit, which requires the energy- 
momentum conservation, we look for the best 
pairing of photons into tt" by constructing a 
pseudo-x^ variable for each of the 15 possible 
pairs: 

^J^l^i ) J = l,2,...,15. 

^=l V ^"VO / 

(4) 

where 

- nij^^o is the invariant mass of the i^^ir^ , in 
correspondence of the j*'* combination; 

- M^o is the 7r° mass, (M^o = 134.98 MeV/c^ 

my, 

- CTm g is the resolution on m^o . 

The chosen pair is that one which minimize the 
X^- A second kinematic fit constraining the 
mass has been performed, improving the resolu- 
tion on 2; by a factor two. The Monte Carlo, MC, 
z distribution at generation (pure phase space), 
see figdJand after reconstruction, see figlH shows 
that resolution effects are not negligible for this 
analysis. Three samples with different efficiency, 
e and purity on pairing, P: 

Low Purity P = 75.4% e = 30.3% 
Medium Purity P = 92.0% e = 13.6% 
High Purity P = 97.6% e = 4.3% 

have been analyzed by cutting on the difference 
of the two lowest value of x^- 

4. Measurement of the slope parameter a 

In order to estimate a an unbinned likelihood 
function is built by convoluting the event density 
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Figure 1. MC z distribution according to pure 
phase space. 



Figure 3. Data-Monte Carlo ratio of the z distri- 
bution. The MC distribution is pure phase space. 




Figure 2. MC z distribution after selection and 
reconstruction. 



simulation we have realized that this is essentially 
due to a different value of invariant mass of three 
pions system (t\P tt" tt" ) in the Monte Carlo gen- 
erator, = 547.30 MeV/c^, with respect to the 
one recently measured [il3 by our experiment: 

Mr, = 547.822±0.005stat±0.069sj;stMeV/c2. (6) 

As a consequence, the accessible phase space on 
data is larger than the one on Monte Carlo simu- 
lation. A new evaluation of a has been obtained 
after a kinematic fit with an additional constraint 
on the 77 mass, see Table [T] 



with the resolution function and correcting for 
the probability of wrong photon pairing in tt*^ 's. 
Using a sample with high purity and fitting in the 
range (0 — 1), we found the preliminary result [T3]: 

a = -0.014 ±0.004 (siai) ± 0.005 (sysi), (5) 

where the systematics has been evaluated by 
varying the analysis cuts, the fit range and mea- 
suring the maximum observed variation of a with 
respect to the samples with different purity. 
We observed a relevant dependence of a on the 
fitting range. Moreover, the low purity sample 
shows two different slopes (see figE]) in the Data- 
MC ratio of the z distribution. With a dedicated 



Range 


Low Pur. 


Med. Pur. 


High Pur. 




a(-10-'^) 


a(-10-3) 


a(-10-3) 


(0,1) 


-30 ±2 


-31 ±3 


-26 ±4 


(0,0.8) 


-26 ±2 


-28 ±3 


-22 ±5 


(0,0.7) 


-26 ±3 


-27 ±4 


-23 ±5 


(0,0.6) 


-30 ±4 


-31 ±4 


-20 ±6 



Table 1 

Fitted results for the slope parameter a for the 
kinematic fit with the 77 mass constrained at 
M,, = 547.822 MeV/c^. 



In this way, good stability with respect to the 
range of fit is observed, and the linearity of the 
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Data-MC ratio of the z distribution has been re- 
covered, (see fig HI . In order to give the final re- 



x'/ndf 21.51/23 
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RES 


-9 


-4 


-3 


Range 


-4 


-4 


-3 + 3 


Purity 


-f + 3 


+4 


-4 


BKG 


0. 


0. 


-I + l 


Mr, 


-f 


-2 


-5 


Total 


-10 +3 


-6+4 


-8+3 



Table 2 

Summary table of systematic uncertainties. The 
total systematic uncertainty is obtained by 
adding in quadrature each contribution. 



Figure 4. Data-Monte Carlo ratio of the z dis- 
tribution after the kinematic fit with 77 mass con- 
straint. 



suit, the phase space area where the z distribu- 
tion is flat, is used as fit region (z e [0, 0.7]). The 
systematic uncertainties on a are summarized in 
Table [21 We have first taken into account the 
correction for the Data-MC discrepancy in the 
photon energy resolution, RES, the effect of the 
fit range. We have then determined the a varia- 
tion with respect to other purity samples and the 
error due to the combinatorial background evalu- 
ation. We have also tested the a-dependence on 
the chosen value of the used in the fit con- 
straint. 

5. Results 

We quote as preliminary result for the slope pa- 
rameter a the one obtained with the Medium Pu- 
rity sample , (about 650000 77 decays). The 
result including the statistical uncertainty from 
the fit and the evaluated systematic error is: 

a = -0.027 ± 0.004 (stat) tomt {syst) (7) 

with: x^/ndf = 13.72/17. The resuh is within 
errors compatible with the result from Crystall 
Ball based on 10^ events of: 



The our measurement of a also agrees with the 
calculations from the chiral unitary approach [5] . 
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